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Abstract We perform n-body simulations for models with a DE component. Besides of 
DE with constant negative w = p/ p > —1, we consider DE due to scalar fields, 
self-interacting through RP or SUGRA potentials. According to our post-linear 
analysis, at z — 0, DM power spectra and halo mass functions do not depend 
on DE nature. This is welcome, as ACDM fits observations. Halo profiles, 
instead, are denser than ACDM. For example, the density at 10 h^^kpc of a DE 
~ W^^ Mq halo exceeds ACDM by ^ 40%. Differences, therefore, are small 
but, however, DE does not ease the problem with cuspy DM profiles. On the 
contrary it could ease the discrepancy between ACDM and strong tensing data 
(Bertelmann 1998, 2002). We study also subhalos and find that, at z = 0, the 
number of satellites coincides in all DE models. At higher z, DE models show 
increasing differences from ACDM and among themselves; this is the obvious 
pattern to distinguish between different DE state equations. 
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1. Introduction 

Deep survey and CBR data confirm that ~ 70 % of the world is Dark En- 
ergy (see, e.g., Efstathiou et al 2002, Percival et al 2002, Spergel et al 2003, 
Tegmark et al 2001, Netterfield et al 2002, Pogosian et al 2003, Kogut et al 
2003), as needed to have the accelerated expansion shown by SNIa data (Riess 
et al 1998, Perlmutter et al 1999). The nature of Dark Energy (DE) is a puzzle. 
ACDM needs a severe fine-tuning of vacuum energy. DE with constant neg- 
ative w = p/ p > —1 has even less physical motivation. Apparently, the only 
viable alternative is dynamical DE, a classical self-interacting scalar field (p 
(Wetterich 1985). Among potentials with a tracker solution, limiting the 
impact of initial conditions, Ratra-Peebles (1988, RP hereafter) and SUGRA 
(Brax & Martin 1999, 2000) potentials bear a particle physics motivation. 
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Figure 1. Spectrum evolution for 

ACDM (solid line) and RP (long dashed); 
the dot-dashed line is the linear prediction 
for ACDM. 
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Figure 2. The virial density contrast Ac 
vs. Q,m at 2 = 0. Models are indicated in 
the frame 



Studying a dynamical DE model requires: (i) a linear treatment, to yield 
CBR spectra and transfer function; (ii) a post-linear treatment, to yield halo 
virial density contrasts and mass functions; (iii) a non-linear treatment. Here 
we report results on (ii) and (iii). We use the n-body program ART, modified to 
deal with any dependence of (matter density parameter) on a (scale factor). 
Mainini et al (2003b) give analytical fitting formulae for such dependence. 
Further details are in Mainini, Maccio' & Bonometto (2003a) and Klypin et al 
(2003). RP and SUGRA are parametrized by the energy scale A/GeV. 

2. Non-linear results 

Fig. 1 shows the evolution of the spectrum, as obtained from simulations 
of ACDM and RP models (A/GeV= 10^), the most distant models treated. 
Models were normalized so to obtain the same number of halos at z = 0. 

Halos were extracted from simulations using the virial density contrasts Ac 
obtained by Mainini et al (2003b), where one can find plots for the dependence 
Ac(a); here we show Ac dependence on at z = (Fig. 2) Figs. 3 & 4 
show the mass function and its evolution in a number of models. 

Using ART facilities, a particular halo was magnified in all simulations. 
Fig. 5 shows that its profile is NFW with a concentration depending on DE 
nature. Concentration can also be considered on a statistical basis. Fig. 6 
shows how halo concentrations depend on the model. Here concentrations are 
defined as the ratio between the radius Tc at which the density contrast is 110 
and the radius in the NFW expression of the radial density. 

We also studied how the number of satellites of a halo depends on DE nature. 
In Fig. 7 we report such dependence. However, also in this case, once care is 
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Figure 3. Mass function at 2 = and 
2 = 2 for the same models of Fig. 1 . Evo- 
lution is faster for ACDM than for RP 



Figure 4. Halo number evolution in var- 
ious models. Unlabeled curves refer to 
SUGRA and constant w = -0.8 




Figure 5. Density profile for a single 
magnified halo. Solid, short dashed, long 
dashed lines refer to ACDM, SUGRA, RP. 



Figure 6. Concentration distribution in 
various models 



payed to properly normalize numbers to the same central halo velocity, no 
appreciable dependence on DE nature can be found. 

3. Conclusions 

In this paper we showed how a simple modification of the program ART 
permits to perform a wide analysis of dynamical DE models. This task is 
simplified by the very structure of the program, which uses the scale factor a 
as time-variable, and requires only 

dt/da = H^^^a VLm{a) /^m{ao) 

{Ho', today's Hubble parameter), to detail the action of forces. Once ilm(a) 
(the dependence of the matter density parameter on the scale factor) is as- 
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Figure 7. Number of halo satellites 

signed, the dynamical problem is then properly defined. Most of the prelimi- 
nary work was then performed at the post-linear level. This provided us suit- 
able expressions for the virial density contrast, so that halos can be selected 
in the correct way in all models, and also the required fitting expressions for 
rim(a). 

Discriminating DE models from ACDM essentially requires good data at 
high redshift. A discrimination at z = can be made only using an observable 
sensitive to the concentration distribution. In principle, such an observable 
exists and is related to strong lensing (giant halo statistic). Further work in this 
direction is in progress. 
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